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1. INTRODUCTION 


THE main features of the structure of collagen at the molecular level have 
been reasonably well worked out (Ramachandran and Kartha, 1955; 1956). 
The polypeptide chains form three non-coaxial helices in the protofibril and 
are further wound up into coiled coils about a common axis. The three 
chains are attached to one another by means of hydrogen bonds approximately 
perpendicular to the fibre axis. These triple chain rods are arranged in a 
hexagonal array, the distance between them varying with the amount of 
moisture present. Though this structure explains fairly well the wide angle 
X-ray pattern of collagen, some of the details are not yet clear. The most 


puzzling aspect is the non-occurrence of the 1120 (110) reflection or its ana- 
logue having the same €-value in the higher layers. It is usually supposed 
that the Fourier transform of the structure is weak in this region, and that the 
absence of a spot in the zero layer corresponding to the 110 reflection is due 
to this cause; but it is difficult to believe that it should be so for every layer 
line. In fact, a spot having roughly this é-value (€ = 0-20) occurs on the 
fourth layer line of kangaroo tail tendon and was identified as such by Rama- 
chandr2n and Ambady (1954). However, the actually observed é-value 
differs from the expected value of 0:24 by an appreciable amount. Thus it 
would appear that the whole series of reflections of the type 11/ is forbidden 
for collagen. Further it is found that the observed é-values in general lie 
close to 0:14, 0-28 and 0-42 in all the layer lines, i.e., in the ratio of 1: 2:3 
and not at intermediate values like 1/3 or +/7 times the basic unit of 0-14. 
This suggests that reflections with mixed indices of the type hk (both h and k 
non-zero) do not occur in the X-ray pattern of collagen. 


Although a few spacings larger than 15 A have been previously reported 
(Corey and Wyckoff, 1936) they had not been confirmed. Photographs taken 
in this laboratory of native cattle achilles tendon showed diffuse maxima at 
about 18 A, 29 A and 46A. These agreed with the observations of Corey 
and Wyckoff. The publication of a beautiful photograph taken by Cowan, 
North and Randall (1955), recording a whole series of sharp reflections in 


Al 363 








364 V. SASISEKHARAN AND G. N. RAMACHANDRAN 


this range brought to light the fact that the earlier observations were done 
at insufficient resolution. Even in this photograph the 110 reflection is miss- 
ing on the equator. The fact that such a large number of sharp reflections 
are observed at small angles shows that the specimen is well oriented and it 
is therefore all the more surprising that no reflections were recorded corres- 
ponding to the é-value of the 110 spacing of the hexagonal lattice. 


This difficulty is resolved and at the same time the occurrence of a large 
number of equatorial long spacings is also explained, if it is assumed that the 
structure of collagen is not built on a regular lattice, but on a cylindrical 
lattice. The idea of cylindrical and spiral lattices is briefly discussed below 
and they are applied to the case of collagen. It is very probable that these 
concepts could be extended to other fibrous proteins as well; some suggestive 
evidence in these cases is also presented. 


2. GENERAL THEORY OF CYLINDRICAL LATTICE 


A number of excellent theoretical papers have appeared on cylindrical 
lattices, their classification and their Fourier transform (Whittaker, 1954; 
1955; Jagodzinski and Kunze, 1954; Waser, 1954). We shall consider here 
briefly the general theory of a simple normal cylindrical lattice. Such a cylin- 
drical lattice is built up of a series of sheets, forming the surfaces of coaxial 
cylinders, the radii of these increasing in arithmetical progression. In any 
small region, it could be imagined as having been derived from a_ regular 
lattice by giving it a cylindrical curvature. Thus if pm is the radius of the m-th 
cylinder and pm the number of lattice points along its circumference, then 
Pm = 49 + ma, dy being the radius of the innermost cylinder and 27pm 
=pmb, where b is the distance between the neighbouring points along the cir- 
cumference. Suppose that the z-axis is along the axis of the coaxial cylinders; 
then any point can be defined by the cylindrical co-ordinates p, ¢ and z. Let 
the points which are nearest to the initial plane (¢ = 0) lying on the m-th 
cylinder have the co-ordinates ¢ = 5m. Then the angular parameters ¢ of 
the points in the m-th cylinder are given by 


= ) 
¢ Pat ™ 


where v is an integer such that 1<v<pm. Also p= pm=d)-+ ma and 
z= nc, where c is the repeat distance along the axis of the cylinder. 


In accordance with the usual theory of diffraction, the amplitude of the 
diffracted wave along a direction given by the reciprocal vector s is the Fourier 
transform of these cylinders, which is given by 
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r= 22 Ze0 rs) 5 


where r is the position vector. The vector s may be called the diffraction 
vector, whose terminus is given by the co-ordinates p*, ¢* andz*. Then, 


r.S = pp* cos (¢* — ) — zz* (2) 
Thus, 


T= J op (5 nez*) a exp = [pmp* cos (¢* — dm,)] (3) 


Replacing the exponential functions occurring in the sum in expression (3) 
by a series of Bessel functions, 


T= 2 exp 4 ncez*) 5 5 iT Jg (= pme*) 


m q=—Co 
x exp [ig (¢* — 3m)] X exp (- 2zia?) (4) 
° Pm 
The last sum, %, is finite only when g = kpm (where k is an integer) and it is 
then equal to Pm. It vanishes otherwise, so that, 


2ni ne 
T= SD exp (= ncz*)\ 3 Pm & Skp, (x pmp*) 


q=—0o 


x exp | ikem (s* +5- im) | (5) 
The first summation, Y, has appreciable values only in the vicinity of planes 


in the reciprocal space, given by z* = / A/c, where / is an integer. These are 
the usual layer lines. In the second summation, », the square of the modulus 


of the functions, averaged over rapid oscillations in 4*, causes all cross 
terms between high order Bessel functions associated with different values 
of k, and cross terms between them and Jo, to vanish. Also, cross terms 
involving different m’s tend to vanish unless k =0. Thus when k 40 
the intensity | T |? is proportional to 

= m > (27 ‘ 

2 ~ Pm dina? ( y Pm ) (6) 


This function has peaks corresponding to p* = k A/b (kK = 1, 2, 3,....), 
thus leading to characteristic bands of large intensity in reciprocal space on 
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cylinders, whose radii are multiples of b*. The intensity of the k-th band 
arises essentially from the Bessel functions of orders + kpm. From the 
expression (6), it is seen that the intensities due to the scattering from the 
component cylinders of the lattice are additive for these bands. These bands 
are known as the ‘ diffuse’ series, since they are accompanied by a tail of 
secondary maxima on the longer p* side. 


When k = 0, however, there is interference between the terms J, (2z/A 
< pmp*) for different pm and the amplitudes due to the component cylinders 
are additive. The intensity |T|* is proportional to 


| 2 Pm Jo (FZ pmp*) (7) 


The intensity exhibits sharp maxima on bands corresponding to p* = h Da 
(A integral). 


2 





Thus, the diffraction pattern can be described as consisting of a series 
of layer lines at € = / A/c, on which two series of reflections occur, one the 
sharp series at £ = h A/a (hintegral) and the other the diffuse series at 

=k A/b. 

The above discussion is confined to the case in which the two-dimensional 
lattice inscribed on the cylindrical layers is primitive and rectangular and has 
one axis oriented perpendicular to the cylinder axis. Whittaker (1955) has 
enumerated and classified the possible types of cylindrical lattices and has 
extended the theory to diffraction by regular cylindrical lattices, with an 
oblique generating lattice and by those belonging to the helical series. A 
particular case of diffraction by a helical structure has also been discussed 
by Jagodzinski and Kunze (1954). These papers may be consulted for a 
more detailed account of the theory. 


3. CYLINDRICAL LATTICE FOR COLLAGEN 


Let us consider a cylindrical lattice built up of protofibrils* all running 
parallel to the axis of the cylinders. Then in any cross-section the appearance 
will be similar to Fig. 1, in which each protofibril is represented by a small 
circle. The distance between successive cylinders is a and if b is the distance 
between the protofibrils on any cylinder, then we have the result 
An = nj — nj, = 2na/b, where nj; is the number of protofibrils on the j-th 
cylindrical surface. This must obviously be an integer and the number was 
chosen to be six for collagen. The choice came from a consideration of the 
symmetry of coupling the cylinders together, and the occurrence of a regular 


* The term protofibril is used for the elementary unit of the structure. This may consist of 
a single polypeptide chain as in the a-helix or a triple chain, as in the collagen structure. 
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hexagonal lattice as an approximation. This immediately imposed a rela- 
tion between a and b, namely, that 27a = 6b and thereby the geometry of 
the structure was completely fixed. 


Fic. 1. Cross-section of a cylindrical lattice. 


However, as one goes towards the centre, there are two possibilities in the 
central region, as indicated by Figs. 2 (a) and 2(b). They are best described 
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Fic. 2. The two possible cylindrical lattice structures of collagen. Fach circle repre- 
sents a protofibril. (@) Structure with one protofibril right at the centre. (5) Structure with 
three protofibrils forming an equilateral triangle at the centre. 


by the relations (i) aj = ja and (ii) aj = (j — 4) a. Thus Fig. 2(a) has a 
protofibril right at the centre, while the innermost core in Fig. 2 (b) consists 
of three protofibrils forming an equilateral triangle. In both cases, the struc- 
ture at the centre corresponds approximately to a regular hexagonal lattice 
and outside, it fits in with the scheme for a cylindrical lattice. The choice 
between the two could not be made on a priori grounds, but was possible 
from a detailed comparison with the equatorial long spacings in the X-ray 
pattern, which indicated that the structure shown in Fig. 2 (bd) is the correct 
one. There is a reasonable fit with the data if the central protofibril is not 
there in the structure shown in Fig. 2 (a), 
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The c-axis of the structure was taken to be parallel to the cylinder axis, 
the repeat spacing being 28-6 A, the value for the protofibril. From the 
properties of the transform of the cylindrical lattice, there will be maxima 
corresponding to = h A/a = ha* and = k A/b = kb* on the equator and 
on all layers defined by § = /A/c = Ic*. In the present case, since a/b = 
= 2m/6 = 1-047, the two series of principal maxima (namely ha* and kb* 
series) practically coincide and only one set of principal reflections occur at 
&=h &,(& =~ A/a). On the other hand, no reflections can occur at values 
of ¢ = ha* + kb*, corresponding to mixed indices of the type h, k, with both 
h and k non-zero. Physically, this can also be explained as being due to the 
fact that, while there is a regular order among the points in any particular 
ring in Fig. 1 and also while the distance between any two circles is a constant, 
there is no correlation between the points in different circles in that figure. 


The points get in step and out of step periodically and so the analogue of a 
hkO plane cannot occur. 


Thus, while reflections of the type h0/ and Ok/ are possible, those of the 
type AkO and Aki cannot occur. In particular, 110 is a forbidden reflection 
for collagen, and this €-value cannot also occur in any of the layer lines. 
In fact, the absence of mixed / and k indices in the diffraction pattern of 


cylindrical lattices has been noticed in the case of chrysotile by Jagodzinski 
and Kunze (1954). 


4. COMPARISON WITH THE OBSERVED EQUATORIAL LONG SPACINGS 


The absence of the 110 reflection and its analogue in the other layers 
thus receives a natural explanation. However, the real evidence for the 
cylindrical lattice is obtained from a comparison of the calculated Fourier 
transform of the cylindrical lattice with the observed X-ray diffraction pattern. 
A number of subsidiary maxima occur at intermediate values of ¢ in between 
the principal maxima. These have been worked out for the structure with 
seven cylindrical sheets [Fig. 2 (b)] and the results are shown in Table I for 
the region from 0 to &). The value of spacings larger than 10 A observed 
on the equator and their relative intensities are given in columns 3 and 4 of 
Table I. The first and second columns contain the calculated spacings and 
intensities. The spacings were calculated taking the most intense reflec- 
tion at 12-6A to correspond to the principal maximum of the first order. 
The number seven was arrived at by comparing the theoretical values with 


observations. The agreement between the calculated and observed values 
is seen to be good. 


The sharpness of the observed maxima indicates that all the fibres in 
the specimen have exactly the same number of sheets. On drying the fibres, 
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TABLE I 


Calculated and observed maxima in the diffraction pattern of a 
cylindrical lattice with 7 sheets 











Calculated Observed 
Intensity 
Spacing (arbitrary Spacing Intensity 
units) 

110 43-6 at 
66 12-3 cf = 
49-5 3-6 49-0 S 
38 3-2 37-0 S 
31-5 1-7 31-0 w 
26-5 2-0) ; 

23-5 1-0) 2A°5 m 
20-5 1-8 20-5 w 
18-5 1-0 18-8 m 
16°8 3-2 17-0 m 
15-4 2-9 14-9 m 
13-9 38-0 13-5, V.S. 
12-6* 93-0 12-6* V.S. 








} 


* Adjusted to be equal. 


this regularity is lost and on rewetting, only a streak joining the centre to the 
principal equatorial spots is observed. A similar streak is also observed 
in other layer lines, notably in the third layer line, showing that there is appre- 
ciable scattering power in the region in between the principal maxima. If 
the Fourier transform of the protofibril does not have a large value at the 
¢-values corresponding to the principal maxima, but only in between, then it 
will record as a diffuse reflection on the particular layer line. The broad peak 
at € = 0-20 in the fourth layer line is probably of this type. 


In order to see the behaviour of the transform beyond the first principal 
maximum, the calculations were continued for the range from that to the 
second principal maximum for the structure with seven cylindrical sheets 
[Fig. 2(b)]. From the cylindrical lattice theory, the intensity distribution 
along the equator, or zero layer, is given by 


1(8) =c| 3 32m — 1) Jo(w (2m — 1) Ela*) |? 


ksi med 


+ 2e[ > z 9 (2m—1)*3*sk (am-1) (3 (2m—1) é/b*) 
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where c is a constant. A graph is shown in Fig. 3 connecting intensity and 
the é-value. The intensity calculations were made at intervals of 7 é/a* 
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Fic. 3. A graph connecting the variation of intensity with ¢ for the structure shown in 
Fig. 2 (5). 
= 3€/b* = 0-05. The intensity variations in the range from 0 to £9 (& ~ A/a) 
is mainly due to the contribution by the first part of the above equation. The 
nature of the second function is such that it has appreciable values only when 
£ approaches k &, and thereafter it is oscillatory in nature. The intensity 
however never falls to zero beyond £9, so that sharp subsidiary maxima would 
not be observed. It must be noted that in general there is an appreciable 
scattering power in the region between any two principal maxima. This is 
a feature peculiar to the cylindrical lattice and does not occur with a regular 
lattice, even if the lattice is finite. With a finite regular lattice, the intensity 
distribution between any two principal maxima is the same and exhibits a 
number of subsidiary maxima, which fall away in intensity on moving away 
from a principal maxima. With a cylindrical lattice, this feature occurs only 
between zero and the first principal maxima. The diffuse intensity in the 
region between the principal maxima increases with increasing é-value. 


There is thus clear evidence that the aggregation of collagen protofibrils 
is based on a cylindrical lattice and that they form cylindrical cryptofibrils of 
diameter approximately 200 A [2x (13-15) x 7].¢ In fact fibrils of this 
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+ The term ‘‘cryptofibril” is used to indicate cylindrical structures of the type discussed, with 
a diameter of the order of 200 4, 
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order of diameter showing characteristic 640 A axial period have been observed 
in the electron microscope. In particular, they are very clearly seen in the 
beautiful photographs reported by Swerdlow and Stromberg (1955) of fibrils 
subjected to mercury at high pressure. These pictures indicate that an ordi- 
nary fibril having a diameter of the order of 500-2000 A is built up of such 
cryptofibrils. Burton ef al. (1955) have found that fine threads, about 
150-250 A in diameter, are produced in the early stages of the degradation 
of collagen by chemical treatment. 


5. THE POSSIBILITY OF ROLL STRUCTURE 


So far, we have discussed the cylindrical lattice structure. However, 
it seems physically more probable that a cylindrically curved layer structure 
would have for its cross-section a spiral, rather than a set of concentric circles. 
In that case, the complete structure can be obtained by rolling up a single 
sheet, in which the protofibrils form a parallel array. 


In any cross-section, the equation to the spiral can be written in the form 
r=ké. As before, the whole structure can be characterised by two dis- 
tances (i) the distance a between the sheets and (ii) the distance b between 
the protofibrils in the sheet. We may assume, to start with, that the sheet 
is rolled about an axis parallel to the protofibrils. If this were not so, the 
protofibrils would take up the configuration of a helix, whose radius is con- 
tinuously increasing. So far, it has not been possible to work out theoreti- 
cally the diffraction pattern of a spiral lattice, even of the simple type. How- 
ever, by an experimental approach described in the next section, it is found 
that the rotated Fourier transform is of the same type as the one obtained 
from a cylindrical lattice. 


The spiral lattice appropriate to collagen is shown in Fig. 4. Here also, 
the relation 27 a = 6b was used for plotting the lattice points. The lattice 


Fic, 4. Cross-section of a roll structure (spiral lattice) appropriate to collagen. 
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points were so plotted that the spiral structure thus obtained broadly 
resembled the one shown in Fig. 2 (bd). 


6. OPTICAL TRANSFORMS OF CYLINDRICAL LATTICE 
STRUCTURE OF COLLAGEN 


The above results of the cylindrical and spiral lattice structures were 
confirmed by means of experiments, made on models with an optical 
diffractometer. The diffraction apparatus used in these experiments was 
the one described by Aravindakshan (1957). As with other optical diffrac- 
tion apparatus, this also does not permit the study of diffraction patterns 
from three-dimensional objects. So investigations were carried out only 
of diffraction from plane circular and spiral lattices. However, for the pre- 
sent purpose, this does not involve any limitations, since the diffraction pattern 
obtained from such a lattice is equivalent to the distribution of intensity on 
the zero layer plane in reciprocal space for the corresponding cylindrical 
lattice. 


The pinhole source used was 40 microns in diameter and the resolution 
was fairly high. Photographic etched masks (Aravindakshan, 1957) were 
used at the beginning of the experiment. Because of the non-uniformity 
of the photographic glass plate, distortions were produced in the image. A 
more serious difficulty was a displacement due to the non-parallelism of the 
two surfaces of the glass plate, as a result of which rotated optical transforms 
could not be recorded. The masks were later prepared using a pantograph 
punch on a scale 0-185 mm./A. Each lattice point was represented by a 
hole 0-45 mm. in diameter. The first minimum of the diffraction pattern 
from a hole of this size occurs at € = 0-75, which is well beyond the range 
(€ < 0-30) of the optical transform useful for study. These holes were 
punched on a blackened photographic film. The drawings, from which the 
masks were prepared, were constructed on a scale nine times larger. In 
order to plot the points of the lattice, their loci (either concentric circles or 
the spiral, as the case may be) were drawn first. The lattice points were then 
set off at equal distances along the curve, their positions being checked and 
corrected at least every 90°. The optical transforms were photographed 
directly on a fine grain photographic plate and enlarged for reproduction. 


A comparison with the X-ray results can be best had when the transform 
is obtained with the mask continuously rotated in the optical diffraction 
apparatus. The transform thus obtained may be called ‘ rotated transform ’ 
and for a satisfactory result, neither accurate centering of the mask nor 
absolute uniformity of motion is required. The only requirement is that 





the 
ac! 











Studies on Collagen—II 373 


the plane of the mask must remain the same during rotation and this is readily 
achieved by using a good ball-bearing. 


(a) Cylindrical lattice—The optical transforms obtained for the two 
possible cylindrical lattice structures of collagen [(Figs. 2(a) and 2(5)] are 
shown in Figs. 5 and 6, along with their rotated transforms. The similarity 
between the two transforms is remarkable. In both cases, the six-fold 
symmetry of the pattern arises from the relation 27a = 6b. The sharpness 
of the maxima, their fine structure and the arch-like distribution of the diffuse 
reflections are all to be expected from theory. 


The two rotated transforms contain rings of large intensity only at 
integral multiples of a unit and not at multiples like 1/3, +/7, etc., which 
would be expected for a simple hexagonal lattice. The distinction between 
the cylindrical lattice and a simple hexagonal lattice of the same dimensions 
is easily seen by comparing their rotated transforms shown in Figs. 6 (b) and 
7(b). The absence of rings corresponding to the reflections 1, 1 and 2, 1 
with the cylindrical lattice is particularly noteworthy and supports the theo- 
retical prediction. 


A number of weak rings (subsidiary maxima) are found in the rotated 
transform of the lattices, and are particularly clear in the central region inside 
the first ring. However, their number and spacings differ, e.g., they are quite 
different for the two types of cylindrical lattice structures discussed here. 
In fact, it was the number and spacings of the subsidiary maxima near the 
central region which enabled the authors to choose the structure shown in 
Fig. 2(b) as the proper one. The spacing of the rings for this structure are 
shown in Table II in the second column. The agreement between this and the 
observed spacings is seen to be good. In the photographs reproduced, only 
the stronger among the subsidiary maxima are clearly seen. 


(b) Spiral lattice —The optical transforms of the spiral lattice appropriate 
to collagen are shown in Fig. 8. The similarity between these and the trans- 
forms of the cylindrical lattices is striking; however, the symmetry is only 
approximately six-fold in this case. In the ordinary transform many of the 
spots making up the pattern do not lie in exactly the same positions on the 
various arches, a fact which clearly arises from the spiral character of the 
lattice. 


As before, the rotated transform contains intense rings only at integral 
multiples of a unit. The spacings of the subsidiary maxima near the central 
region obtained from this lattice are given in Table II in the third column. 
The agreement between this and the observed spacings is again good, A 
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TABLE II 


Observed spacings (in A) in the X-ray diffraction pattern of collagen 
and the spacings calculated from measurements on the diffraction 
pattern of cylindrical and spiral lattices in an optical diffractometer 





From rotated transform of 
From X-ray pattern 








Cylindrical Spiral 
lattice lattice 
101 99 
- 67 65 
49-0 49-0 49 
37-0 38-0 39 
31-0 32-0 32 
v (27-0 27 
24:5 124-0 23 
20-5 21-0 20 
18-8 “4 a 
17-0 17-0 17-3 
14-9 15-4 14-8 
13-5 14-2 14-0 
12-6 12-6 12-6 





comparison of the rotated transforms of the cylindrical and spiral lattices 
shows that there is more of intensity in between the principal maxima with 
the spiral lattice. This is to be expected, since there is less order in a spiral 
lattice than in a cylindrical lattice. The phenomenon is seen particularly well 
in the region between the first and second principal maxima in Fig. 8 (d). 
Thus all the evidences available so far are in accord with spiral lattice also, 
with the added advantage that the formation of a roll in the form of a con- 
tinuous sheet is much more reasonable. 


It is quite possible now that collagen has definitely been shown to have 
a cylindrical lattice structure, that other fibrous proteins also may have such 
a structure. This is particularly to be expected since the number of spots 
found in their X-ray pattern are very small, e.g., in keratin, compared to those 
found in the patterns of some polypeptides. In fact, the published X-ray 
pattern of sea-gull feather keratin (Bear and Rugo, 1951) clearly indicates 
the occurrence of a cylindrical lattice. Only row lines corresponding to 
h , are found with h = 1, 2, 3, 4 and €) = 0-045. The absence of any other 
intermediate row lines of sharp spots clearly shows that a cylindrical lattice 
structure is present and also that the ratio a/b is probably the same as in 
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collagen. There is also an extension of reflecting power in some layer lines 
in between the principal row lines. Probably the a-keratin group also has 
this structure, as only one type of repeat spacing has been found along the 
equator. 

SUMMARY 


It is suggested that the structure of collagen is based on a cylindrical 
lattice. The protofibrils, consisting of triple chains, are packed together in 
a hexagonal array close to the centre, and are extended outwards in the form 
of cylindrical sheets. Each sheet has a pseudohexagonal symmetry about 
the common axis. This cylindrical lattice explains (a) the large number of 
equatorial reflections of long spacings observed in the diffraction pattern of 
native fibres and (b) the absence of the 110 reflection and its analogue in all 
layer lines. The number of sheets in a single cylindrical rod is about 7 and 
its diameter about 200 A. 


The above results were verified by model optical diffraction experiments. 
Also, from these experiments it is shown that the structure can also belong to 
the roll type with the added advantage that the formation of a roll in the form 
of a continuous sheet is more probable from the point of view of fibrogenesis. 
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V. SASISEKHARAN AND G. N. RAMACHANDRAN 


EXPLANATION OF PLATE X 


Transforms of the structure shown in Fig. 2a. (a) Ordinary transform, (6) Rotated 
transform. 


Transforms of the structure shown in Fig.26. (a) Ordinary transform, (5) Rotated 
transform. 
Transforms of a simple hexagonal lattice. (a) Ordinary transform, (6) Rotated transform. 


Transforms of the spiral lattice shown in Fig. 4. (a) Ordinary transform, (5) Rotated 
transform. 
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§ 1. INTRODUCTION 


Tue matrix A = (ay), i,k = 1, 2, 3,...... of real or complex constants 
transforms any given sequence x = {xp} of real or complex numbers into a 
sequence y = {A;(x)} by the equations 


(T) yi = Ai (x) =2 dikXk tS & & Opa ). 


The sequence x is said to be summable by the matrix A to the value / if all 
the sums y; are finite and lim yj = /. We denote the transformation (T) 


ico 
by the matrix equation y = Ax, where y and x are column vectors whose 
elements in the i-th row are y; and x; respectively. For the sake of conveni- 
ence we shall use the same symbol to represent a sequence and the column 
vector Corresponding to it. 


Let (co), (c) and (m) denote as usual the spaces of null sequences, conver- 
gent sequences and bounded sequences respectively, with the usual norm. 
If (X) and (Y) are two sets of sequences then (X) (Y) will denote their inter- 
section, i.e., the set of sequences (points) common to both (X) and (Y); and 
we shall denote by J’ (X, Y) the class of all matrices A which have the property 
that xe(X) implies Axe (Y). We define further, for every complex number 
p that I"(c, pc) is the class of matrices A such that lim x;=/ implies 


ico 


lim J ajxxX~ = pl. Therefore in particular I’ (c, c) is the set of all K-matrices 
ico k=1 


(conservative or convergence-preserving matrices) and I'(c, lc) is the set 
of all T-matrices (permanent or regular matrices). 


The summability field of A, i.e., the set of all sequences x = {xn} such 
that Axe(c) will be denoted by the symbol (A), and if (X) is a given set of 
sequences, we define R,x) (A) to be the “* range of A over (X) ”, i.e., as the set 
of sequences y such that y = Ax where xe (X). 
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The object of this paper is to extend some results that Banach [(1), pp. 90- 
95] has proved for I’(c, 1c) to other classes of matrices and to derive therefrom 
a number of other theorems. The immediate extensions of Banach’s results 
and the extensions obtained therefrom of some results due respectively to 
Hanai (3), Hill (4) and Ramanujan (7) are given in §2; in §3 we show that 
applications of these functional analysis methods and result yield a number 
of theorems due respectively to Copping (2), Wilansky (8) and the author 
(5), (6). The proofs as well as some of the results of § 3 are believed to be 
also new. 

§2 

We begin by recalling that the necessary and sufficient conditions in order 
that the matrix A should belong to (i) I'(m, m), (ii) I'(c, c), (iii) T’ (co, cp), 
(iv) I'(c, pc) are, respectively, 


(i) (@): |All shud. F | ain | <0 
1 <i<oco k=1 
(ii) (a) is satisfied, 


(b): lim > dik = p exists, and 


i-co k=1 


(c): lim aj, = 5% exists for each k = 1, 2, 3,.... 


i—co 
(iii) (a) and (c) are true, with 5, = 0 for all k = 1, 2, 3,.... 
(iv) (a), (6) and (c) are true, with 8, = 0 for all k = 1, 2, 3,.... 
THEOREM 1.1!.—If the matrix Ae«I’(c,c) satisfies the condition 
x (A) =lim J ain — J lim ain = p — J bn <0 (2.1) 
i~>co n=1 n=1i~-co n=1 
and if y = {yj} is a convergent sequence such that 
2 ayj=0 (2.2) 
whenever 


E |aj| < co and F ajain =0 for alln = 1, 2, 3,.... (2-3) 


i=1 i=1 


1 Banach [(1), p. 91] has proved this theorem for the case AeJ” (c, 1c) when (2.1) is automati- 
cally satifised. 
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then there exists for every number « > 0 aconvergent sequence x such that 
| Ay(x) — yg | < € for all i=1, 2, 3,.... (2.4) 
ProoF.—The conclusion of the theorem means that ye Rie) (A), where 


Ric) (A) is the closure of the set of A-transforms of convergent sequences. 
By a lemma of Banach [(1), p. 57] it is plainly enough to prove that every 
linear continuous functional? vanishing over Ric) (A) will vanish at y. Now, 
any linear continuous functional defined over (c) is of the form [(1), p. 66]. 


f(9) = coo + E cnn; || +E |en| <0 (2.5) 


where s = {sn} has s for limit. Therefore any linear continuous functional 
over Rie) (A) will be of the form 


F(v) =f (Aw) = Co [pte + 2 8p (up—ue)] + E en E anpitp (2.6) 
where u = {uj} and lim uj= wu». Using the condition (2.1) and (ii) (a), 
we may write ns 

f(v) = F Upbp + Colo x (A) 
where 


bp = Coop + E cndnp for all p= 1, 2, 3,.... (2.7) 


Now if the functional f vanishes over Ric) (A), then taking u!™ = {3;"}, 
where 5,” has the value unity if i = n and is zero otherwise, for n = 1, 2, 3 
....Successively, we see that 


> 


by = 0 for all p = 1, 2, 3,.... (2.8) 
Therefore form (2.6) it follows that f(v) = cots x (A) = 0. 
Since by hypothesis x(A) #0, we see by taking uw = {5,9 that c, = 0. 
The equations (2.7) and (2.8) now give 


¥ cndnp = 0 for all p=1, 2, 3,.... 


n=1 


* We prefer this terminology to Banach’s “linear functional” in order to avoid any possible 
confusion. 
A2 
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and since by (2.5) Y | cn | < co it follows from the hypothesis that Scan = 
n=l n=1 


and therefore f(y) = colim yj + zr Cnyn = 0. This proves the theorem. 
i‘—~>co n=1 
A trivial modification of the proof of Theorem 1.1 gives the following: 
THEOREM | .2.—Jf AeI" (co, Cg) and ye(cy) are such that (2.3): E | a| 
i=1 


< ©, F ain = 0 for all n = 1, 2,.... implies (2.2): ¥ ay; = 0, then there 
i=1 


i=1 
exists for every number « > 0 a null sequence x such that | Aj(x) — yi| < 
Rv of t= |, 2, 350-6 


THEOREM 2.—(Generalised converse of Theorem 1.1.) 


Let AeI’(m, m) and y = {yj} be such that for every « > 0 there exists a 
sequence xe(m) satisfying the condition 


| A(x) — yi | < € for all i=1, 2, 3, 
Then for any sequence {a;} such that 


© \aj| < 00 and ¥ ajain = 0 for all n=1, 2, 3,.... (2.3) 


we shall have a 
z ajyi = 0. 
ProoFr.—Let {a;} be a sequence satisfying the condition (2.3). Then 
a (s) =F ays; is a linear continuous functional defined over (m) and || a || = 


i=1 


y \a;\. Also, ye(m) since Ax = {A;(x)}¢(m) and therefore we have: 


i=1 


; co co 
|a(Ax — y) | = | 3 a4(2 ainxn — ya| 
f=1 n=1 
<|lalj-|Ax—yll<llalj- e (2.9) 
co co , P 
and since Y J |a|-|a@in| > |xn|<llell- |All: xl] < co 
i=l n=1 


we have 


a(Ax) = J a J ainxn 
i=1 n=1 


co co 


= 2 Xn J in 


n=1 i=1 


=0 by (2.3). 


3 Hanai (3) has proved this theorem for the case AeJ"(c, Ic). 





The proof given here is simpler. 
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Therefore (2.9) gives 


|S ays |< ja (Ax — y)| + la(Ax)| <|lal|-«. Since «>0 
4=1 


co 
is as small as we please, we see that Y a;y; = 0. 


i=1 


DEFINITION.—If the matrix A is such that for every sequence {a;} the 
condition (2-3) implies that a; = 0 for all i = 1, 2, 3,...., then A is said to 
be of type M. 


THEOREM 3:1.4A—Let AceI’(c, c) be of type M and let x(A) #0. Then 
Rie) (A) is dense in (c). 

The proof is immediate after Theorem |.1. 

As a consequence of Theorems |.2 and 2 we have the 

THEOREM 3-°2.°—The matrix AcI'(Co, Co) is of type M if and only if 
Ric, (A) is dense in (co). 

THEOREM 3.3.—Let AeI’(m,m). Then Rie) (A) is dense in (c) if and 


only if the set of sequences (points) consisting of the columns of the matrix A 


co 
and the sequence of row-sums {3 ain} is fundamental in (c). 


n=1 


Proor.—If we denote by A the linear space generated by the set of se- 
quences {3;'}, n = 0, 1, 2,....., where 5; = 1 or 0 according as i=n 


or i <n for n> 0 and 8; = | for all i, then it is well known that A is (c). 
Now, the linear space generated by the set mentioned in the theorem is the 


same as the set Ra (A). Since A is (c) and since A defines a linear continuous 
operation over the Banach space (m) [in virtue of (1), Theoreme I, p. 54]and the 


relations /.u.b. | dinXn|<l.u.b. J | ain | + lu.b.|xn| = || All: || x} 


l<i< co n=1 l<i<co n=1 l<n<co 


it follows that Ra (A) > Rye) (A). Therefore we have: 








Ra (A) = Ra (A) > Ree (A). 


* Banach [(1), p. 93] has proved this for AeJ”(c, 1c) when the condition x (A) 4 0 is auto- 
matically satisfied. 

5 It may be remarked that Ramanujan (7) has proved this theorem for AeJ” (c, 1c) without 
however making use of the condition (ii) (+) of §2. 
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But R~ (A) C Rye) (A) and hence we see that Ra (A) = Ryo (A). This 
proves the theorem. 


By arguing in exactly the same manner we get the result: Let A«I"(m, m), 
Then Rye, (A) is dense in (co) if and only if the columns of A form a set 
fundamental in (co). Theorem 3.2 now yields the following: 


THEOREM 3.4.—The matrix AeI"(co, Co) is of type M if and only if its 
columns form a set fundamental in (co). 


THEOREM 4.—(Generalised converse of Theorem 3.1.) 
If AeI’(m, wm) and Rim (A) is dense in (co) then A is of type M. 
ProoF.—Since Rim) (A) > (co) we see that for every ye (co) and every 
number « > 0 there exists a bounded sequence x such that 
| Ay (x) — yx | < € for all i= 1, 2, 3,.... 
In virtue of Theorem 2, the condition (2.3) will then imply that Y ay; = 0. 


Since this is true for every null sequence y, we see by considering the sequences 
y™ = {5}, n = 1, 2, 3,....that aj = 0 for all i= 1, 2, 3,.... and the 
theorem is proved. 


COROLLARY.°—The matrix AeI(c, c), x(A) #90, is of type M if and 
only if Rye) (A) is dense in (c). 


THEOREM 5.1.—If AcI” (co, Co), then for every bounded sequence % 
such that Axge(Co) and for every number « > 0 there exists a null sequence 
x satisfying the condition 


| Ag (x) — Aj(%) | < € for all i=1, 2, 3,.... 


In other words, if AeI" (co, Co) then Rye, (A) is dense in Ria.) «m (A), 
where (Ag) denotes the set of sequences which are summable by A to the value 0. 


Proor.—Let us write y; = Aj(X) for each i= 1, 2, 3,.... and let 


{a;} be a sequence such that > | aj | < coand 2 a44in = 0 for all n = 1, 2,3, 


1=1 t=1 


We have then 





6 Hill (4) has proved this for reversible matrices in J"(c, 1c) though he does not seem to use 
the reversibility. 





(th 


R, 


sul 


all 


nu 


ye 
th: 











Applications of Banach Functional Methods to Summability 383 


co co co co 
DS aye = FT Aji(X%) = LY 4 LT AinSn 
i=1 i=1 i=1 n=1 


co co 
= J Sn = in = 0 
i=1 


n=1 


(the order of summation being justified by the fact that 2 ¥ | ag || din | 
i=1 n=1 


‘[snl< Z| a4] + All + Isl] < 0), 


The result follows now from Theorem 1.2. 


THEOREM 5.2.°—IJf AeI"(c,c) and y(A) 40 then Ree) (A) is dense in 
Ryaycm (A). 


Proor.—Same as for Theorem 5.1, only at the end we use Theorem 
1.1 instead of Theorem 1.2. 


THEOREM 5.3.—Let AelI’(m,m) and let there exist for every number 

« > 0 and for every given xye(A)(m) a sequence x such that 
| Ay (x) — Aq (%) | < € for all i=1, 2, 3,.... Then for any 
sequence satisfying (2.3) we shall have ¥ aA; (xo) = 9. 
i=1 
Proor.—We have only to write y = {yj} = {Aj (x»)} and apply Theorem 2. 
§3 

In this Section we apply the results obtained above to reciprocals of 

summability matrices. 


THEOREM 6.1.—Let AeI(c,c), x(A) #0 and let 'AeI’(m,m) where 
‘A.A =I, the identity. Then (A)(m) = (c). 


PRooFr.—We note that ‘A sums every xeRi) (A) for ‘A (Au) = u for 
all we(c), indeed for all we(m). Now, by Theorem 5.2, Rie) (A) is dense 
in Ryay¢m (A); also, since ‘AeI’(m, m) we see that ’A defines a linear conti- 
nuous operation over (m) by a well-known result [Banach (1), Theoréme 1, 
p. 54]. Therefore it follows that ‘A sums every xeR,a)ym)(A). That is, if 
ye(A)(m) and x = Ay then ‘Ax = ‘A(Ay) = y belongs to (c). It follows 
that (A) (m) c (c) and since Ael"(c, c) the theorem is proved. 

THEOREM 6.2.°—/f AelI'(c, c) and x(A) =0, then there is no matrix 
‘AcI’(m, m) such that 'A.A =I. 


7 Banach [(1), p. 93] has proved this for the special case AeJ"(c, 1 c), 
$This result has recently been announced by Wilansky (8), 
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PROOF.—It is well known that under the hypothesis of the theorem 
there exists a bounded divergent sequence summable by the method A, i.e, 
(A) (m) #(c). The required result now follows from Theorem 6.1. 


THEOREM 6.3.°—If AcI"(c, c) and 'AeI"(m, m), 'A.A =I, then (A) (m) 
== (c). 


ProoF.—This follows immediately from Theorems 6.1 and 6.2. 


THEOREM 6.4.°—Jf AeI’(c, c) and A 'eI’*(m, m) is such that A.A 
= A.A!=I, then A'clI'(c, c). 


ProoF.—Let x = {xn} be any convergent sequence. Then A (A-!x) = 
xe(c). Theorem 6.3 gives now that A-'xve(c). Since this is true for all 
xe(c), the theorem is proved. 


THEOREM 6.5.—If AcI” (Co, Co) and ‘Ae’ (m, m), ‘A.A =I, then 
(Ao) (m) = (Cp). 


The proof is the same as that for Theorem 6.1, only read (co) wherever 
(c) occurs and use Theorem 5.1 instead of Theorem 5.2. 


THEOREM 6.6.—/f AeI (co, Co) and A“«I’(m,m) is such that A.A 
= A.A‘ =I, then A-*el"(Co, Co). 


The proof follows from Theorem 6.5 just as Theorem 6.4 followed 
from Theorem 6.3. 
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CERTAIN results regarding the asymptotic values of sums involving the 
fractional parts of numbers have been proved and made use of by PILLAI 
and CHOWLA® in some other connection. Here we shall prove three theorems 
about sums involving the fractional parts of numbers, but of a different type. 
Throughout this paper {x} denotes the fractional part of x and [x] denotes 
the integral part of x, where x is any positive real number. 


THEOREM 1: J 3 = (1 — y) x + 0(x*), where y is Euler’s constant. 


d= 


THEOREM 2: z d {*} = (1 = *) x? + 0(x log x). 


THEOREM 3: zd a Fe } = (.- {kt ) xk 4 O(x*), for all 


k > 2, where €(s) is the Riemann Zeta Function. 


PROOF OF THEOREM 1: We have 
E d(n) =xlogx + (2y — 1)x +0(x), (1) 


(vide HARDY and WRIGHT,! Theorem 320), where d(n) is the number of 
positive divisors of n including | and n. From the definition of d(n) we have 


d(n) = 
Hence 

2 d(n) = D> oe 

n=1 din 
= 1 which is the same as 

dd’<x 

e @ [eld 

Pd@ =F fF 

nz d=1 d’=1 


-Z lal 
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= ECG {a) 


=X (1og x+y+0 (<)) — > {3}, (2) 


d=1 


since 


z t=logx ty + (1). 


x 
From (1) and (2) we have the theorem. 


PROOF OF THEOREM 2: We have 


2 oc (n) = as a 0 (x log x), (3) 


(vide HARDY and WRIGHT,! Theorem 324). Here o(n) is defined as the sum 
of the positive divisors of n including | and n, and so 


a(n) = zz d. 
dia 
Hence 


Fom=F 70a 


n=l n=1 din 


= » d 
dd’<x 

z (eid) 
es ee 


d=a a’=1 
= 24(5- {3}) 

=x.x— 3d ]3}. (4) 
So from (3) and (4) Theorem 2 follows. 


To prove Theorem 3 we shall require the following lemmas: 


LemMA 1: (The Euler-Mclaurin sum formula in its simplest form, 
vide KNopp,? p. 529.) 


If f (x) is continuous and differentiable for all values of x, then 


Sfx) = f f(dx+4(SO +f) 


+ f (x — [x] — Pf’ (x) de. 





) 








Applying the above lemma to the case when 


i 
f(x) _— (x + 1)e 5: CE 0, we have 
LEMMA 2: 
n l ni-@ ; - _ 
2 x _ l—a oa O(n ), as 0. 
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PROOF OF THEOREM 3: As usual, let us denote by ox, (nm) the sum of the 
k-th powers of the positive divisors of m including | and n, for all k >2. 


So we have 


ox (n) = Sd. 


din 


We shall now estimate S' ox, (n) in two different ways. Firstly, 


n=1 


ZS on(n) =F 5 dk 


n=1 n=1 din 
= J dk 
dd’<x 
{zid’] 


=> dk x I 


d=1 @/=1 
=2(0- tas) 


2 1 @ 
=x 2 ga- Z dt i 


d=. d=1 


Since 1 — k <0, using Lemma 2 in the right side of the above, we have 


2 k 
xy ox (n) =x G. oo 0 (x*-1) — 
= ioe + 0(x*) _ x 
d=1 
Next we have 
ZS on(n)= 5 dk 
n=1 dd’<x 
a2 {eida’] . 
= Piz d+ Fz! 
d’=1 dz1 d’=@/2 


= § + S’ (say). 


“t. 


d=1 


“f 


(eid) 
= a* 


qd=1 


PERE: 


. 


(5) 
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Since — k <0, by Lemma 2 we have 


[x/d’}+k 


S= F (Fo toed’). 
But 

ET =(-of-0<«< 

“GY 0-9) 

and in 

1 <d’ <5 
. < 4 and so 

ed’ < 4 


x 


uniformly and so 


le] =G) tr +0()} 
ELE GY 
“PBato(eBa) 


@’=1 d’s1 


Hence 


Since k > 1, using the expression for the Riemann Zeta Function for the 
sums in the above, we have 


k+1 
S= 64 (EK +1) +00) + 00x49) (2) + 004) 


_o(kK+1) x k 
a es x*+1 4. Q (x*), 


F [eid’} 
Now in x|2 <d'< x, x\d’<2 and hence in the sum ¥ d*, there are atmost 
d=1 


[eid’) 
two terms. So we have + d*=0(1) for any a, when x/2<d'<x 


dei 
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Hence 


s'=0( 3 


== §(»x). 
Hence 


Fon) =S +S’ 


n=l 


= EME ED xk + O(x4) + O(n) 


a C(k + 1) K+1 _ k 
= ay Mt + OC. (6) 
Now equating the right sides of (5) and (6), we get the theorem. 


Finally, I wish to thank Professor V. Ganapathi Iyer for his valuable 
advice in the preparation of this paper. 
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ABSTRACT 


Radiations emitted in the decay of Samarium-153 have been studied in 
the Siegbahn-Slitis beta-ray spectrometer. Using the internal conversion 
electron spectrum and the photo-electron spectrum with tin as radiator, 
the internal conversion coefficient a, has been determined for 102 Kev. and 
70 Kev. gamma-rays. The relative intensities of the three beta-ray branches 
have been determined. A weak gamma-ray of energy 83 Kev. has been found 
and can be interpreted as a transition to the ground state from the first 
rotational level in Europium-153. The multipole order and character 
of the 102 Kev. transition is discussed. 


INTRODUCTION 


SAMARIUM-153 (half-life 47 hrs.) decays by B- emission to several levels in 
Europium-153 and has been studied by various workers.'-® There has been, 
however, considerable disagreement about the internal conversion coefficient 
aj,2>4.9.10 for the 102 Kev. gamma-ray and some uncertainty about the relative 
intensities of the various beta-ray groups involved in the decay scheme. 
As a result of experiments on Coulomb excitation, rotational levels at 82 Kev. 
and 187 Kev. have been established in Europium-153 and it was of some 
interest to find out whether these levels are also formed in the beta-decay of 
Samarium-153. 


Accurate determination of internal conversion coefficients has assumed 
special importance as recently there has been a suggestion,!*-!* with some 
experimental evidence™ in its favour, that in the case of magnetic dipole 
(M,) transitions, the coefficient has a smaller value than that predicted theo- 
retically by Rose,!® on account of the finite size of the nucleus. The 102 Kev. 
transition in Europium-153 has been considered by previous workers!:2:4 
to be a mixture of M, and E, types and the coefficient a; as determined in 
this work is discussed from this standpoint. 

The method of investigation adopted in this work consisted of three 
main experiments. In the first place, the continuous beta-spectrum was 
recorded and the Fermi-Kurie plot was constructed. Secondly, the internal 
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conversion electron spectrum was investigated, giving the relative intensities 
of conversion lines for different transitions, K/L ratio and the Auger lines. 
Finally, using tin as radiator, the photo-electron spectrum was _ recorded, 
giving relative intensities of unconverted part of the transitions and X-rays 
in Europium, consequent on internal conversion. From these measurements, 
the conversion coefficient was calculated. This method of using photo- 
electron spectrum has not been employed by previous workers for this iso- 
tope and has the further advantage of detecting gamma-rays of weaker inten- 
sity and of low energy. 
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EXPERIMENTAL OBSERVATIONS AND RESULTS 


A source of Samarium-153 of high specific activity was obtained from 
A.E.R.E., Harwell. Fig. 1 shows the continuous beta-ray spectrum of 
Samarium-153, N/Hp being plotted against Hp, the momentum. Fig. 2 
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is the Fermi plot, the continuous curve having been obtained from the 
observed points. Assuming allowed shape for different beta-groups, this 
is seen to be made up of three befa-ray groups with maximum energies at 
792 + 10 Kev., 685 + 15 Kev. and 610 + 20 Kev. respectively. The rela- 
tive intensities of the beta-groups and the Log ft. values for these transitions 
are given in Table I. 














TABLE I 
Group Energy in Kev. % Transition log ft. 
By 792+10 20 7°3 
Be 685+15 65 6°5 
Bs 610-+20 15 71 





Fig. 3 shows the internal conversion electron spectrum. The source 
was deposited on a thin aluminium foil (0-13 mg./cm.”), and was about 
100 wg./cm.” in thickness. Five electron peaks are clearly observed and are 
interpreted as arising from two transitions corresponding to 102 Kev. and 
70 Kev. and Europium X-rays, as shown in Table II. Absorption in the 
counter window, which was 80 yg./cm.? thick, was taken into account, using 
correction factors obtained from auxiliary experiments using various values 
for window thickness. 

















TABLE Il 
Line Sum Gamma 

energy Interpretation Relative intensity (Kev.) energy 

(Kev.) (Kev.) 
20°3+-5 .. K,(Kconv. of 70 Kev.) 0-29 +- -03 69-5 70 
33+1 * Ay (Auger, K-2L) AytAg sanere ks 
38+2 .. Ag (Auger, K-L-M) K, 193: 08 
53-341 .. K,(Kconv. 102 Kev.) I 102-5 102 
62 .. L,(L conv. 70 Kev.) 0-06-+- -01 70 70 





L, (L conv. 102 Kev.) 0-16+ -02 101-5 102 
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Photo-electron spectrum of gamma-rays emitted in the decay of Sama- 
rium-153 was studied in the present investigation, using different radiators, 
namely, tin, uranium and gold. Fig. 4 shows the photo-electron spectrum, 
for tin as radiator. In this experiment, the source was kept close to an alu- 
minium disc, thick enough to absorb all primary beta-rays. A layer of tin, 
250 wg./cm.? was deposited by evaporation, on one side of the disc, facing 
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the spectrometer entrance slit. The line energies measured from the 
spectrum and their interpretation are shown in Table III. Lines due to 
gamma-rays at 102 Kev. and 83 Kev. are observed, besides the lines due to 
Europium K, and Kg X-rays. The line due to 70 Kev. gamma-ray, which 
is highly converted internally, coincides in energy with the L photo-line of 
Kg X-ray which appears as a weak line on the high energy wing of the 
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intense L photo-line due to K, X-ray in Europium. The relative intensities 
of the 70 Kev. and 83 Kev. gamma-rays, in comparison with 102 Kev. gamma- 
ray were found to be -06+ -01 and < -04 respectively, from the photo- 
electron spectra as well as the scintillation spectrum. 


TABLE III 
(Radiator: Tin) 





Intensity ratio 























Line energy Interpretation Sum N,/N, after 
(Kev.) (Kev.) necessary 
corrections 
10-5+-8 .. Kconv. Eu K X-rays 39-7 0:77+0-05 
16°5+1 .. Keconv. EuK X-rays 45-7 
19+1.. .. Auger line in Tin K-L-L 
35-5+1 .. Leonv. Eu K X-rays 40 
5542 .. .. Konv. 83 Kev. gamma 84-2 
73-5+1 .. Kconv. 102 Kev. gamma 102-5 
9742 .. .. Leonv. 102 Kev. gamma 101-5 





©ig.6 shows the decay scheme of Samarium-153. Three beta-groups 
leading respectively to the 170 Kev. 102 Kev and ground states in Europium 
and the two main gamma-transitions at 102 Kev. and 70 Kev. are shown. The 
weak radiations at 170 Kev. and 530 + 20 Kev. which have been reported 
by previous workers*-*»® and at 600 Kev.® were observed neither in the 
photo-electron spectrum nor in the internal conversion spectrum clearly and 
therefore only an upper limit can be put for their intensities as 1%. 


However, they were observed in the scintillation spectrometer with 
relative intensities < 1% in comparison with 102 Kev. gamma-ray. Since 
Samarium oxide (Sm,Og) ‘ spec-pure’ used for neutron irradiation contains 
normally Europium as a small impurity, there is a need for some further 
investigation on these weak transitions. However, contribution due to 
these weak transitions to the X-ray intensity will be negligible on account 
of their small internal conversion coefficient, and it has been neglected in 
the present case. The weak gamma-ray at 84 Kev. is shown to be due to the 
de-excitation of the first rotational level above the ground state. Europium 
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155 (half-life 1-7 years) has a gamma-ray at 87 Kev. and can be produced 
as an impurity in the sample but its activity will be less than a:thousandth 
of that of Samarium-153. Also the relative intensities of 83 Kev. and 102 
Kev. radiations are seen to remain unchanged over a period of ten days, 
indicating that 83 Kev. gamma-ray. observed here is due to Samarium-153. 
For the purpose of estimating a,, the conversion coefficient of 102 Kev. 
gamma-ray, and for estimating the relative intensities of the beta-branches, 
the 83 Kev. transition will not be.considered. It is also assumed that all 


the decay is by B- emission and that there is no decay by electron-capture. 
A3 
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In that case, the vacancies in the K-shell of Europium, giving rise to K, wh 
and Kg X-rays, will all be due to internal conversion of gamma-rays y, and De 
Yq Of energies 102 Kev. and 70 Kev. respectively. 
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Fic. 4. Photo-electron Spectrum 
Nz, the number of K-shell vacancies is given by 


Ng = Nek + Neke 


where Nex, and Nex, represent the number of gamma-rays y, and yz con- 
verted in the K-shell. 


Or, 





Nx = Nek, (1 + ng 


Nex, (1) 


Now, 
Nz = N,or 
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where Ny is the number of X-rays and w, is the fluorescence yield. 
Denoting the number of unconverted transitions by N,, 


Nx = Nr. @k 
Ny, Ny, 
N (1 + a) 
Nx —_ Nxt = _— i . Nek,/ 
Ny,%k N,, Ny, 
By definition, a,,, the internal conversion coefficient of y, is 7 
: 
N 
- 1 ») 3 
cae ax, ( * Nek, ( ) 


The value of ax, is given by equation (3). The term Nex,/Nex, on the right 
is the ratio of intensities of K-conversion lines due to y, and y, and is 
obtained from the conversion electron spectrum (Fig. 3) Nz/N,,, on the left 
is the ratio of the intensity of X-rays and of unconverted y,-rays and is 
obtained from the photo-electron spectrum, from the area of the photo-line 
due to y, and the area under all the photo-lines due to X-rays (Fig. 4), 
correcting for the variation of photo-electric absorption cross-section in 
tin with energy. Correction has also to be applied for absorption in the 
aluminium absorber used for absorbing primary beta-rays. One obtains 


Nz _ Npg @#y,° (1 — eT y,!) 

Ny, Noy, e-Mgt qd — et, Tal) 
where Npz and Np,, are the areas under the photo-peaks due to X-rays and 
yrTays, jy, and yz are total absorption coefficients for y,-rays and X-rays 
in aluminium, thickness ¢ (440 mg./cm.*) and wz, and wry, are corresponding 
photo-electric absorption coefficients for tin, thickness, (0-250 mg./cm.2), 
I, as obtained from the tables.” 


The fluorescence yield w, is obtained from the internal conversion 
spectrum, from the equation 


Nz (1 — Wk) — Nae. 


Where Nx, the number of K-shell vacancies is given by the area under the 
K-conversion lines of the transitions corresponding to y; and y, and Nae, 
the number of Auger electrons, ~ the area under the Auger lines. The 
value of w_ is found to be 0-87 * ° ya in agreement with wy, = -88 + -02 
given by K. Siegbahn.? a, for 102 Kev. gamma-ray as calculated in the 
above manner is found to be 0:69 + -08 which agrees with -66 + -1 obtained 
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in the usual manner from the internal conversion line and the f spectrum 
on which it is superposed. These present values can be compared with 
ay, = 0°62 + -15 given by Lee and Katz* and a,~-6 given by Siegbahn.? 
Now a, for 70 Kev. transition can be calculated by using the a, for 102 
Kev. transition, the ratio of K internal conversion, i.e., Nex,/Nek, and the 
ratio of gamma-ray intensity for these two transitions; a, for 70 Key. 
transition comes out to be 3-5+ 1-4. This value of a, for 70 Kev. can 


be compared with 3-8 given by McGowan,” and 4-4 + -4 given by Dubey 
et. al.® 








(+x) THE INTERNAL CONVERSION COEFFICIENT 
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Fic. 5. Photon Energy in Units of mc.? 


The branching ratio for the f, transition, leading to the ground state 
in Eu®, is readily found from the Fermi plot (Fig. 2) and is 20%. One can 
calculate the branching ratio for B, and f; transitions, making use of the 
internal conversion data found from the conversion electron spectrum. 
Referring to the decay scheme and denoting the number of y, and yz transi- 
tions by N, and N, respectively, and the number of the two beta-transitions 
by Ngs and Ng,, one gets, 


Ne, = N, 
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and 
N, = Na, + Na, 
Ni Net Np 1 4 Na : 
Ng Ng, Ne, ‘ 


Also, the number of transitions N, equals the number of unconverted 7, 
rays, N,,, and the number converted in K and L shells. 


N, 7 Nek, + Nex, + N,, 


Ns = Nex, (1 + yt + yt 
L 1 
Ni = Nex, (1 + + 5 
similarly, 
L 1 
N: = Nex, (1 + K, + ian 


Li 1 
Ni _ New E+ + ae, 


N; Nex (,,l2, 1)" 
e+e < 





The quantities on the right, conversion coefficients, K/L ratio and the ratio 
of intensities of the K-conversion lines for y, and y, are known from the 
internal conversion spectrum. Hence, N,/N, and Ng,/Ng, can be calculated. 
The branching ratios for B, and B; are found to be Ngo/Ng, = 5. 


The results are summarised below: 














Gamma-ray a K/L Intensity 
102 Kev. 0-69+0-08 6-2+ -8 1 
+1-3 
70 Kev. 3-5+1-°4 5— 9 0-06+0-01 
83 Kev. > ae < 04 





By/Be/Bs = 9°3/1/0-2 
DISCUSSION OF RESULTS 


One of the objectives of this work was to determine the K-conversion 
coefficient for 102 Kev. gamma-ray, about which there has been a lack of 
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agreement among previous workers. McGowan!® and Marty® found 
it to be 1:14+ -2 and 1-2 +0:1 respectively, while Lee and Katz* and 
K. Siegbahn? reported it to be of the order of 0-6. The value as determined 
here is 0:69 + -08. In the light of this, it is of interest to consider the 
character and multipole order of the 102 Kev. transition. Previous workers 
have assumed this transition to be a mixture of E, and M,. Table IV below 
gives the theoretically expected values for conversion coefficient a,, and K/L 
ratio, for this energy and character E,, E,, M, and My. Fig.5 gives the 
value of the coefficient a, for Z = 64, against energy in units of mc.? down 
to -3, for E,, E,, E, and M, and M, types of transition and the curves as 
extrapolated to 0:2 corresponding to 102 Kev. These curves are plotted 
from the tables of conversion coefficient by Rose ef al.,!*> and the values for 
a; in Table IV are obtained from these extrapolated curves. 


TABLE IV 





E Z2jE a, K/L 





E, E, M, M, E, E, M, M, 
102 Kev. 38:53 0:25 0-9 15 90 60 20 $793 5-0 








It will be seen from the above consideration that the experimental 
value for the conversion coefficient for 102 Kev. transition, i.e., 0-69 + -08 
is not in agreement with any one of the theoretical values for the four types 
of transitions. The value, however, would be consistent with a mixture 
of E, and M, types with a preponderant E, component. 
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The other alternative is that the transition is a mixture of M, and E, 
types, with a correction factor operating in the case of the M, coefficient, 
due to the finite nuclear size, as recently discussed by Sliv!? and Church.}* 
In that case, one finds that it is a predominantly M, type of transition and 
that the correction factor for the coefficient is of the order of 0-5, somewhat 
larger than that to be expected from Sliv’s calculation for nuclei in this 
range of atomic number. 


From shell model considerations the assignment (d;,. +), (ds). + ) and 
(sy. +) respectively to the ground, 102 Kev. and 170 Kev. *® states seems 
resaonable and this would favour the second alternative (M, + E,) for the 
102 Kev. transition rather than (E, + M,), as the latter requires a change 
of parity. 
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§ 1. INTRODUCTION 


AMONG the most beautiful and interesting phenomena in crystal optics are 
the interference figures displayed by crystalline plates in convergent light. 
Figures of the most diverse descriptions are met with even in the realm of 
transparent crystals, and this is not a matter for surprise: for, firstly the 
crystal may possess optical activity—in which case the two waves propagated 
along any direction in the medium are in general elliptically polarised; 
and secondly the state of polarisation of the incident light as well as that of 
the light for which the analysing arrangement is completely transparent may 
be adjusted to be linear, circular or elliptic in form. By the use of the 
Poincaré sphere, the varied phenomena which fall under this head may be 
brought under a general method of analysis, which—without complicating 
the explanation of the simpler phenomena—helps in the understanding of 
more complex phenomena. 


As specific applications of the method we shall consider some of the 
interference figures displayed by the optically active crystal, quartz. Using 
a single basal section, an interference figure consisting of two spirals coiled 
around each other may be obtained with a circular polariser and a linear 
analyser. This case has been theoretically discussed in the treatises of 
Walker,’ Pockels? and Szivessy* (though, due to certain minor errors, the 
final description of the spirals given in these treatises differ as regards the 
handedness and orientation of the spirals). When a usual polarising micro- 
scope is used, the quarter-wave plate has be to inserted just prior to the 
analysing nicol, and hence the spirals have to be observed under an arrange- 
ment slightly different from that considered in the above-mentioned references. 
The spirals thus observed with a linear polariser and a ‘ circular analyser ’ 
are discussed in the treatise of Mascart.t The sense of description of the 
spirals is reversed when the quarter-wave plate is transferred from a posi- 
tion behind the quartz plate to a position in front of it (though Mascart 
states that the same results are obtained under both arrangements). 
402 
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We shall also discuss from a novel point of view, the Airy’s spirals 
observed when two basal sections of quartz of equal thickness and of 
opposite sign are superposed, and the combination observed between cros- 
sed nicols in convergent light. 


§2. DISCUSSION OF THE GENERAL CASE 


As is well known, each point in a convergent light figure is the focal 
point of a particular bundle of parallel rays emerging from the plate, of 
which each pair of rays which have been derived by the splitting of a common 
incident ray, can (after passing through the analysing arrangement) 
interfere at the focal point. Hence our first step is to discuss the interference 
effects presented in parallel light at normal incidence by an arbitrarily cut 
plate of any transparent material when examined between an elliptic pola- 
riser C, and an elliptic analyser C,.f The crystal being transparent, the 
two waves propagated through it must necessarily be in two oppositely 
polarised states A and A’ as indicated in Fig. 1 (since only then will the 
intensity of the emerging beam be independant of the phase retardation intro- 
duced by the plate [Part I,° §2)). 


When the incident light (of unit intensity) in the state of polarisation 
C, is split into two beams of opposite polarisation A and A’, these component 
beams will have the intensities cos?4a,’ and sin? 4a,’ respectively, where 
a,’ is the angular distance between the states C, and A on the Poincaré 
sphere (Part I, §2). The analyser C, transmits the fractions cos* 4a,’ and 
sin? 4a,’ of these beams, since it transmits only the corresponding resolved 
components®»® of these beams in the state of polarisation C,; here a,’ is the 
angular separation between the states C, and A on the Poincaré sphere. 
Taking the state of relative phase of the two oppositely polarised beams at 
the point of entry in the medium as the zero or standard, their phase dif- 
ference at the point of emergence from the crystal will be equal to 5 (which 
is the path retardation suffered by the slower beam of polarisation A’ rela- 
tive to the other). The phase difference between the corresponding resolved 
components transmitted by the analyser will however be smaller than this 
value by an angle A; according to Part I, eq. 9, A denotes now half the area 
of the lune AC,A’C,A, or A = ZC,AC, = ZC,A‘C,. (The angles are 
measured with the usual sign convention, being reckoned positive in a 
counter-clockwise sense [see Fig. 1]). 





+ An elliptic analyser C, appears completely transparent for light of elliptic polarisation C,; 
a quarter-wave plate followed by a linear analyser together form a simple elliptic analyser (Part I 


§8). 
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Thus the analyser C, transmits two beams having intensities I, and I, 
and phase difference A given by 


I, = cos* 4a,’ cos? $a,'; I, = sin? 4a,’ sin? 4a,’ \ 


These beams being in the same state of polarisation, the final intensity I 
transmitted by the analyser will be given by the usual interference formula 


I=1,+1, + 2VII, cos A (2) 
The expression for the intensity I can also be written as the difference of 


. . . . on 
two terms, the first of which is the intensity cos? 4 C,C, which would be trans- 
mitted if the plate were absent (5 = 0), the second term giving rise to sub- 
traction colours in white light. 


Equations (1) and (2) are the basic equations on which the further dis- 
cussion of all interference phenomena shown by transparent crystals are 
based. Though they have been derived for the most general case they are 
practically identical with those customarily used’ for the simplest case, viz., 
when the plate possesses only linear birefringence, and the polarising and 
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analysing elements are nicols—in which case A is equal to zero or 7. In 
discussing any specific case, the values of the angles a,’, a, and A which 
obtain in that case have to be used; but their explicit values need be 
introduced not in the general expression (2), but only in certain subsidiary 
expressions determining those particular aspects of the phenomena which 
are to be investigated. 


The variation of I with the direction of propagation may be observed 
in convergent light, using a plate cut normal to an optic axis—each point 
in the convergent light figure corresponding to a definite direction of pro- 
pagation. For a plate of moderate thickness the retardation 5 introduced 
by the plate increases rapidly as we proceed outwards along directions 
normal to the curves of constant retardation. The rate of variation of 
I as we proceed outwards from the optic axis is therefore taken to be 
predominantly due to the change in A. We should therefore expect the 
appearance of interference rings, the curves of minimum intensity occur- 
ring along directions for which the pairs of pencils transmitted by the 
analyser ‘ destructively ” interfere, i.e., along the curves where A = (2n + 1) 2. 
These do not however coincide with the curves where the retardation of the 
plate is an odd multiple of 7. The physical reason for this is that the final 
phase difference A between the interfering pencils transmitted by the analyser 
is not identical with the retardation introduced by the plate but exceeds it 
by an amount (— A); this latter may be correctly described as the contri- 
bution to the phase difference arising from the processes of dissolution and 
analysation, and it depends on the mutual relation between the state of 
incident polarisation C,, the states of polarisation A and A’ of the beams into 
which it is split, and the analysing state C,. The curves of minimum inten- 
sity are thus given by 


S§=(2n+1)r+A (3) 


and from what has been said above, these curves of minimum intensity do 
not in general follow the curves of constant retardation, 5 = const., because 
the angle A is not constant but is itself a function of the direction of pro- 
pagation (since the states A and A’ vary with the direction of propagation). 
A fairly accurate method of plotting the curves of minimum intensity (for 
a uniaxial crystal or a biaxial crystal of not too small an axial angle) is 
obviously the following; the circular curve 6 = (2n + 1)7 is first drawn 
and the radii vectors of this curve are increased by amounts which vary with 
the azimuth and which correspond to an additional retardation of A. In 
this manner it can be shown that the interference rings exhibited by a basal 
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section of quartz between inclined nicols take the form of squares with 
rounded corners. 


The intensity at any point on a curve of minimum intensity may be 
obtained by substituting A = 7 in (2) and is given by 


Tain = COS* $ (4y’ + Gy’) (4) 


This is not constant along a curve of minimum intensity, since the states 
A and A’ vary with direction, and hence the rings appear darkest along 
particular zones which have to be determined for each specific case. 


§3. SpmaL FiGures DUE To A SINGLE BASAL SECTION OF QUARTZ 


In discussing the convergent light figures of quartz we will have to keep 
in mind the following features regarding the propagation of light near the 
optic axial direction. Referring to Fig.2.a, consider any point P in the 
convergent light figure, whose polar co-ordinates are r, 9 (the origin O re- 
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presenting the optic axial direction, and OV being a vertical line). Along 
the direction represented by P are propagated two crossed elliptic vibra- 
tions; and since quartz is a positive crystal, the major axis of the slower 
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vibration coincides with the plane of the radius vector OP itself (so that 
@ is the azimuth of the major axis of the faster vibration A with respect to 
the horizontal). Hence if ¢ be the ellipticity of the faster wave, its state of 
polarisation will be represented on the Poincaré sphere (Fig.2 5), by the 
point A of longitude 20 and latitude 2«, « being positive or negative accord- 
ing as the crystal is optically left or right handed. The numerical value of 
the ellipticity decreases rapidly as the angular distance r from the optic axis 
increases, (tan 2¢ = 2p/Ar*), being a maximum along the optic axial 
direction where opposite circular vibrations are propagated. In what 
follows the main fact to remember is that for all points P on a circle of radius 
r, the ellipticity is constant—so that the latitude 2¢ of the point A is constant, 
only the longitude 26 being altered. 


We shall first discuss the interference figure presented by a left-rotating 
basal section when a left-circular analyser (represented by the pole C; in 
Fig. 2 5) is used, the vibration direction of the incident plane polarised light 
being vertical (represented by the point V). In this case we have A= ZC,AV. 
It is seen immediately from Fig. 2 b that (keeping 2¢ constant) as 20 increases 
from 0 —> 40 —>7 — 32/2—> 27, the ZCjAV also increases continuously and 
goes through the same range of values 0-+4n—-27; —2—->—4r->0. 
This means that along the circle 5 = (2n + 1) 7, the phase difference (A bet- 
ween the interfering pencils transmitted by the analyser falls short of 
(2n + 1) (condition for destructive interference) by an angle C,AV which 
increases continuously with the azimuth 9 — and which in fact becomes 
exactly equal to twice the azimuth at = m7/4. It follows immediately 
(see Section 2) that the curves of minimum intensity consist of two mutually 
enwrapping left-handed spirals which are related to one another by a rotation 
of 180°. From the triangle VAC; we have 


tan ZC,AV = tan 26/sin 2¢ (5) 


In the usual treatment* equations (3) and (5) are derived by a more lengthy 
procedure, and these equations may be used to discuss the form of the 
spirals in more detail. Close to the optic axial direction we have 
LC\AV ~ 20, so that the spirals are given by: 


8 = V/(Ar? + 4p?) = (2n + I) 7 + 20 


If the two arms of the spiral be extrapolated to the origin (where they 
actually fade away) the common tangent at the origin will be at an azimuth 
(47 + p) with respect to the vertical—where p is the total optical rotation 
along the optic axis. Towards the border the figure must approximate to 
the non-spiral form shown by inactive crystals. The transition occurs by 
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way of the non-uniform rate of increase of the arm of the spiral—which 
manifests itself as kinks, which in turn assume the proportion of disconti- 
nuities towards the border of the figure (see Pockels, Plate VI, Figs. 3 and 4). 


The zones along which the successive arcs of the spiral appear darkest 
may be determined with sufficient accuracy by studying the variation of 
(4) along a circle described about the optic axis (a, = constant). Referring 
to Fig. 2b, the arc a,’ acquires its maximum value of (7 — 2) at 20 =0, 
and its minimum value of 2¢ at 20= 7. Along directions close to the optic 
axis (where 2e > 7/4) the sum (a,’ + a,’) lies between the limits 47 to z. 
Hence the arcs of the spiral appear darkest along the vertical radius of the 
field of view. This result has also been derived in the usual presentation. 
We must emphasize however, that for a plate of moderate thickness (say 
3mm.) the above result holds only for points on the first convolution of 
the spiral, since at greater angular distances 2« < 7/4 (see Pockels, Plate VI, 
Figs. 3 and 4). 


When the analysing state C, is changed to its opposite state C,’, the 
result (as is seen from Fig. 1) is to diminish the value of A by 7 and to change 
the length of the arc a,’ to its supplement. The same alteration in the ex- 
pression (2) may be produced if the polarising state C, is changed to its 
orthogonal state. Hence, when a right-circular analyser is used, the entire 
spiral figure is rotated around by a right angle compared to the previous 
case (when a left-circular analyser was used). 


We may obviously summarise the above results in a form applicable 
to both right and left-rotating basal sections. The handedness of the 
double spirals exhibited with a linear polariser and circular analyser is 
always the same as the handedness of the quartz. When the handedness 
of the circular analyser is opposite to that of the quartz, the tangent to the 
spiral at the origin coincides with the vibration direction of the light 
emerging along the optic axial direction; and, close to the optic axis, the 
spirals appear darkest along the diameter perpendicular to the plane of 
vibration of the incident light. A change in the handedness of the circular 
analyser causes a rotation of the entire figure, as such, through a right angle. 


We see from Fig. 1 that when the state of polarisation of the incident 
light, and the state of polarisation C, of the light for which the analyser is 
transparent are interchanged, the sign of A is reversed, but the remaining 
factors in (1) remain unchanged. In Fig. 2 b the same result is obtained by 
change of the sign of 26. Hence the entire spiral figure exhibited with 
a circular polariser and linear analyser is merely the reflection about the plane 
of vibration of the analyser, of the figure obtained when the polariser and 
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analyser are interchanged. This statement holds with regard to all the 
details of the figure, and hence we need not discuss this case further. 


The results of this section were qualitatively verified using a plate of 
right-quartz. The general arguments may be applied mutatis mutandis 
for discussing the spiral figures in negative uniaxial crystals and in biaxial 
crystals. In the latter case only a single spiral is formed because, as the polar 
angle of the point P increases, the azimuth of the major axis of the faster 
elliptic vibration increases at only half the rate. 


§4. Arry’s Sprrats DUE TO Two SUPERPOSED BASAL SECTIONS 


Our general treatment can also be applied to the case of superposed 
transparent plates which often exhibit curious interference figures (see 
e.g., Walker, Joc. cit., p.293). For this purpose it is necessary to replace 
the effect of the passage through a succession of plates by a passage through 
a single ‘equivalent plate’. This can be done by combining successive 
rotations of the Poincaré sphere. 


As an example, let us suppose that a left-handed and a right-handed 
section of quartz (both of the same thickness and cut at the same angle to 
the optic axis) are superposed (in that order) such that the corresponding 
principal planes of the two plates are in coincidence, the combination being 
then viewed in parallel light at normal incidence. 


Referring to Fig. 3, let the point A, of latitude 2€ represent the state of 
polarisation of the faster elliptic vibration (of ellipticity €) propagated in the 
first plate; the state of the faster elliptic vibration of ellipticity (— ¢) pro- 
pagated in the second plate is then represented by a point A, of the same 
longitude as A,, but of latitude (— 2¢). Let us construct an isosceles trian- 
gle A, X A, as indicated, such that the base angles A,A,X and XA,A, are 
both equal to 5/2 where 4 is the retardation of each plate. 


The action of the first plate is equivalent to an anti-clockwise rotation 
of the sphere about A, through twice the internal angle at A,, while the 
action of the second plate is equivalent to an anti-clockwise rotation about 
A, through twice the internal angle at A, (see e.g., Part I, § 6). By a well- 
known theorem for compounding rotations,® the two successive rotations 
are equivalent to a single rotation about the equatorial diameter through 
X, through twice the external angle at X. We have thus proved that the 
combination is equivalent to a single optically inactive birefringent plate of 
retardation A, the faster vibration direction of which makes an angle — a 
with the common principal plane of the two quartz plates (which corrtes- 











Fic. 3. 


ponds to the major axis of the faster ellipse). Here a and A have the values 
indicated in Fig.3. From the spherical triangle AyXX,) we have 


tan 2a = tan $6 sin 2 (6) 


Since the equivalent plate does not possess optical activity we should 
expect the appearance of ‘ isogyres’ when two equal basal sections of left- 
and right-quartz are superposed, and the combination observed between 
crossed nicols. The ‘isogyres”’ should occur along the directions where 
the ‘equivalent’ principal planes of the combination coincide with the 
vibration directions of the polariser and analyser. Since, as we have seen, 
the equivalent principal planes for any direction of propagation do not 
coincide with the principal planes of the individual plates the ‘isogyres’ will 
not take the form of a uniaxial cross. If 6 be the azimuth of a point in the 
field of view with respect to the plane of the polariser, the dark ‘isogyres’ 
obviously occur where @ is equal to a or (47 + a), where a is given by (6). 
This equation has been derived by an entirely different method of analysis 
in the usual presentation (Walker,’ p. 368, eq.39). According to this 
equation it may be shown by following the usual treatment that the dark 
‘ isogyres’ take the form of four mutually enwrapping left-handed spirals 
known as Airy’s spirals. Besides these spirals we will have dark curves, 
where the retardation A of the equivalent plate is a whole multiple of 2z, 
and from Fig.3 these occur along the circles § = 2mm. 


The sense of description of Airy’s spirals. is reversed when the right- 
handed plate is placed first because the sign of a is then changed. 
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5. SUMMARY 


The geometric method of specifying states of polarisation by points 
on the Poincaré sphere is used for giving a unified and physically intelligible 
approach to the interference phenomena displayed by crystalline plates in 
parallel or convergent light—under general conditions when the polarising 
and analysing states are linear, circular or elliptic in form. Examples dis- 
cussed are the spiral figures exhibited in convergent light (a) by a_ basal 
section of quartz between a circular polariser and linear analyser (or vice- 
yersa), and (b) by two superposed basal sections of left- and right-quartz 
between crossed nicols. The Airy’s spirals observed in the latter case are 
interpreted as the ‘isogyres’ of the optically inactive plate to which the 
combination is equivalent. 
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ABSTRACT 


In investigations preliminary to the study of the electrolysis of certain 
rare earth metal and thorium salts in anhydrous basic solvents, it was 
found that anhydrous acetates of Lat’, Cet+® and Th** are moderately 
soluble in anhydrous ethylenediamine. The resulting solutions show the 
conductance behaviour of weak electrolytes. Cerium (III) acetate solution 
in ethylenediamine undergoes oxidation in air to cerium (IV). 


INTRODUCTION 


THE electrodeposition of the rare earth metals and thorium is usually accom- 
plished by the use of their molten halide mixtures; but the possibilities of 
lower temperature depositions from electrolytic solutions are interesting from 
the theoretical and practical aspects. Inasmuch as the highly electropositive 
nature of these metals preclude their electrodeposition from aqueous solutions 
of its salts owing to the preferential hydrogen discharge, such procedures 
must of necessity involve basic non-aqueous solvents. Such solvents, because 
of their strong donor characters, provide media in which solvated electrons 
can readily exist and in which lower oxidation states are strongly stabilised. 
The stabilities of alkali and alkaline earth metals in liquid ammonia are well 
known. The formation of solvolysis products of limited solubility with 
anhydrous rare earth metal and thorium compounds preclude the use of 
liquid ammonia. The favourable electrolytic properties of anhydrous ethy- 
lenediamine,? the success which has been achieved in the electrodeposition 
of certain alkali? and rare earth metals from this medium** and the failure 
to electrodeposit thorium from thorium (IV) chloride in this solvent! prompted 
investigation of electrolytic behaviour of anhydrous acetates of some rare 
earth metals and thorium in anhydrous ethylenediamine. The present investi- 
gation reports data on solubilities and conductances of acetates of lantha- 
nuin (IIT), cerium (IID) and thorium (IV) in anhydrous ethylenediamine. In a 
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subsequent communication, results of electrolysis experiments will be de- 
scribed. 
EXPERIMENTAL 


1. Purification of solvent—Commercial ethylenediamine monohydrate 
(B.D.H.) was dehydrated and purified as previously outlined! The an- 
hydrous product had a specific conductance of 1-35 x 10-* mho cm.“ at 30° C. 
It was stored in glass-stoppered flasks sealed with Dow-Corning silicone high 
vacuum grease which has been found to have very little effect on electrical 
conductance in this solvent. These flasks were kept in desiccators contain- 
ing soda lime, magnesium perchlorate and phosphorus (V) oxide and were 
opened only in a dry box since ethylenediamine absorbs both moisture and 
carbon dioxide from the atmosphere. 


2. Preparation of metal acetates—The usual method of obtaining an- 
hydrous metal acetates was by the high temperature reaction of the respective 
oxides with excess quantities of ammonium acetate, followed by the removal 
of the unreacted ammonium salt by heating first in a current of dry nitrogen 
and then in vacuo.® Anhydrous acetates of lanthanum (IID), cerium (IID 
and thorium (IV) were obtained by the reaction of the respective analytically 
pure nitrates with acetic anhydride.® It is found that this method adopted 
by Panda and Patnaik® is easy and results in purer products than that obtained 
by Moeller and Zimmermann through the above high temperature reaction.® 
The composition of these anhydrous acetates were established by ignition to 
oxides in tared platinum crucibles over a Meker burner. 


3. Apparatus and manipulation—All anhydrous materials were handled 
in moisture- and carbon dioxide-free air in a “Gloved Box, Berkeley 
Type,” Model—RI1104 supplied by Scientific Service Inc., California. All 
the solutions were prepared in this dry box and all apparatus containing such 
solutions was sealed by Dow-Corning silicone high vacuum grease before 
removal to the atmosphere. 


4. Solubility determinations—The respective powdered metal acetate 
was added in small amounts slowly and with brisk shaking to 20 ml. of an- 
hydrous ethylenediamine in a big test-tube fitted with ground glass stopper 
until some sizable amount remained undissolved. Portions of this suspension 
kept in this stoppered test-tube were thermostated for five days at 30+ 0-05° C. 
with frequent agitation. Then the tubes were centrifuged for twenty minutes 
and opened inside the dry box. The rare earth and thorium contents of these 
solutions were determined by precipitating the hydrous hydroxides from 
measured volumes with excess water and weighing as ignited oxides in a plati- 
num crucible. 

AS 
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5. Conductance measurements.—Conductances of the various solutions 
were obtained after thermostating at 30 + 0-05°C. at definite dilutions as 
previously outlined,’ using a Philips Conductivity Bridge, Model GM4249/01, 
the current switch of the Bridge being set on 1000 cycles per second and a 
conductance cell fitted with ground glass stopper and bright platinum elec- 
trodes. The cell constant was 0:2022cm.-!_ The specific conductance of 
the solvent was subtracted from that of the solution in each case. 


RESULTS AND DISCUSSION 


Solubility determinations—Solubility data obtained are summarised in 
Table I along with the solubility data obtained by Moeller and Zimmermann® 
for purposes of comparative study. 


TABLE I 


Solubilities of some rare earth metal and thorium salts in 
ethylenediamine at 30° C. 






































Solubility in | 
Salt ia Ref. | 
g./100 ml. moles/litre g./100 g. | 
of soln. of soln. of solvent | 
La(C,H,O;), .. 6-97 | 0-22 | 8-01 | | 
La(NO;); .. si “A | 2-45 | 9 
Ce(C,H,0,); .. 10-64 0-34 | 12-41 
Nd Br; m J a | 0-94 9 
Nd I; E ” " 2:53 9 
Nd (C,H;0,)5 .. 7-81 0-30 9 
Y (C,H,0,); .. 6-26 0-30 9 
Th (C,H,0,), .. 3-73 0-08 4-30 | 
Th Cl, = 11-08 0-30 = = 
Hg (C,H,0.), | 0-19 ul 








Anhydrous lanthanum (IID) acetate and thorium (IV) acetate gave colour- 
less solutions in ethylenediamine while the solution of cerium (III) acetate was 
deep brown. That this discolouration was due to atmospheric oxidation of 
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some part of cerium (III) to cerium (IV) was revealed during the dilution of 
this solution with water when an immediate violetish blue precipitate of 
hydrous hydroxide of cerium (IID) and cerium (IV) was obtained. Similar 
atmospheric oxidation of cobalt (II) salts to cobalt (III) salts in ammoniacal 
solution is common in co-ordination chemistry. The moderately high solu- 
bilities of these anhydrous acetates are probably due to the relatively higher 
covalent nature when compared with the corresponding anhydrous metal 
halides and nitrates as apparent from Table I. The enhanced solubility of 
thorium (IV) chloride may be the result of extensive solvation.! Of the an- 
hydrous rare earth metal acetates so far investigated, cerium (IIT) acetate has 
the highest solubility in anhydrous ethylenediamine. 


Conductance studies—Typical conductance data (equivalent conduc- 
tance, A mhos cm.? equiv.—! and concentration of the salt in gram equivalents 
per litre, C) as summarised in Table II and Fig. 1, indicate that the an- 
hydrous acetates of lanthanum (III), cerium (III) and thorium (IV) behave as 
weak electrolytes in anhydrous ethylenediamine. 


TABLE II 


Conductance data for ethylenediamine solutions at 30° C. 


























| Lanthanum (III) acetate | Cerium (III) acetate Thorium (IV) acetate 
Komen | vero] x (Con | vent] x flamem S| ve atat| a 

| 0+3308 575 2-614 | 1-0060 1003 1-816 0-1040 | 323 1-215 
| 0-1654 407 2-896 | 0-5030 709 1-886 0-0520 | 228 1-435 
| 0-0827 288 3+306 | 0+2520 , 502 1-978 | 0-0260 161 1-814 
| 0-0414 203; «3-941 | 0-1260 355 2-271 | 0-0130 114 2-312 
| 0-0207 144 | 4-720 | 0-0630 251 2-594 | 0-0065 81 3-076 
| 0-0104 102 | 5-526 | 0-0315 178 3-156 | 0-0033 57 3°748 
| 0-0052 72 | 6-668 | 0-0158 126 3-917 | 0-0016 40 4-791 
| 0-0026 51 | 8-021 | 0-0079 89 4-974 0-0008 29 6-073 
| 0-0013 36 9-883 | 0-0039 63 5-962 0-0004 20 9-414 
| 0-0007 26 | 11-240 | 0-0020 43 8-078 . ae oe 




















Among the three salts investigated, the strength of the electrolyte in 
ethylenediamine increases in the following order ;—thorium (IV) acetate, 
cerium (III) acetate and lanthanum (III) acetate. The observed order 
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of increase in equivalent conductance parallels the increase of their 
ionic radii (Th** — 1-10; Cet? — 1-18; Lat®— 1-22°A) with the con- 
sequent increase in their ionic nature. A comparison of conductance of 
anhydrous thorium (IV) chloride! and anhydrous thorium (IV) acetate shows 
that the former is a slightly stronger electrolyte in ethylenediamine than the 
latter since the latter salt is obviously more covalent than the former. The 
conductance behaviour of other anhydrous salts of the trivalent rare earth 
metals*.5.*. and a number of other salts!® indicates that they behave as weak 
electrolytes in ethylenediamine. This is not surprising when one remembers 
that compounds containing highly charged cations dissolved in a solvent of 


only moderately large dielectric constant (ethylenediamine at 25°C. is 12-8) 
behave as weak electrolytes. 
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§0. INTRODUCTION 
THE equations of motion of a viscous liquid, viz., 
> 
4 — q x curlg + }grad q? = — grad (7 +2)+ dg (0.1) 
involve non-linear terms and several devices (inverse or semi-inverse methods) 


have been introduced for studying the equations (0.1). The notion of 
additivity or superposability of two motions is one such device that has 


. . tad 
received considerable attention in recent years. If (q,, p,;, 2,) are the flow 
velocity, the mean pressure and the external force potential of a viscous 


+> 
liquid motion and (q2, pz, 2.) are the similar quantities for another motion 
of the liquid, the two motions are additive resulting in the motion 


C4 + qos Py + Pe +7, 2,, + 2,) if and only if?” 


+ > + a+ 
curl (q; < curl gg + qe X curl q,) = 0. (0.2) 


The adjusted pressure 7 is determined from the equation 
> > > > TT > > 
4, X curl gg + q. + curl q, = grad ¢ +h: q). (0.3) 


The addition of plane motions of viscous liquids has been examined 
in detail by Ram Ballabh* and the problem for axially symmetric motions 
with the flow velocity vector in the meridian plane has been studied by 
P. L. Bhatnagar and P. D. Verma.** In the present paper we examine the 
problem of addition of axially symmetric motions for which the flow 
velocity vector does not necessarily lie in the meridian plane, so that the 


oe . 
velocity g can be represented in the form 
+ _ ~> 


> 
q = Qxix + In! + Gule (0.4) 


~ 
w 


*I am grateful to Prof. P. L. Bhatnagar for bringing this to my notice, 


418 








moti 
besic 
addi 
repr 


lines 
clas: 


to | 
mot 
the 
circ 


iS | 


an 


T 





oe 


ds) 
of 
has 


OW 
us 


ion 
on 





Addition of Axially Symmetric Motions of Viscous Liquids 419 


in the notation of L. M. Milne-Thomson.t The axial symmetry of the 
motion requires that gy, In Io depend only on the co-ordinates x and & 


besides the time ¢. We characterize all such steady motions which are self- 
additive (Theorem 2 below) and see that the flow velocity can be explicitly 
represented in a special type of such flows. 


REMARK |. The class of (axially symmetric) motions with the stream- 
lines coincident with the vortex lines is included as a special case of the 
class of (axially symmetric) self-additive motions. 

REMARK 2. The condition for a steady motion of a viscous liquid 
to have the circulation-preserving property is the same as that for the 
motion to be self-additive. The results of § 2 below thus also characterize 
the class of all axially symmetric, steady, viscous liquid flows, which are 
circulation-preserving. 


§1. Conditions for addition of motions in the case of axial symmetry. 
if 


+. 
is the velocity vector of a motion, the continuity equation, viz., divg =0 
is satisfied on writing 


_ 1 8% _ 1 by 
fe & 8a? Ta ~ & 8x 


and we may write q, = x. 


Taking the two motions given by 


ae = (a oe) e+ (OSE) iy tre HD (LD) 


Ww 


we obtain from condition (0.3) that 


8 8 i 
3 Ae [se EM) 4. so (E* *b) + + 5 xx) | Ix 


1 


81 (ps 8 
+=, [Sat (Es) + 5 a = m (E%,) + = (@ *xuxd | ip, 


1 dy, 5 ~ dy. 5 dy, 8 ~ 
+ or? [se 5x (@ X2) +> dar 5x (& Xv pe 5x Soy (& x) 
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~ HF ig @r] | 
= grad (= +41 - ax) (1.2) 


where E? is the differential operation 


std 8 
a2 3a? & 8a 


> = 
Writing 7/p + q, - dg = F, we have the right-hand side of (1.2) in the form 


=. o> 1 3F > 
Sx le + 5 in +e Sq a (1.3) 


> > 
Since q,, q2 are axially symmetric fields and 7 also is to be axially symmetric, 
we have iat = 0. In view of this, we get from (1.2), (1.3) 


ae | ge EW) + SEM) + = @ a, dx 


+ oe [ Se CE) + So EMA) + gre @? xixd | dex 


dF 


F dx e+ de = dF at (1.4) 
and 
5 (41, D Xe) 5 (ve, & X1) ae 
3(x,a) > 80,0) ~° (1.5) 


By means of cross-differentiation with respect to x, we may derive from 
(1.4), the single equation 


5 (Ys, E’ bs 2) | 5 (Yo, E*y;) 
8 (x, a) . 5 (x, &) 


— 2 [EM + SE) +2 @uxd] =0 1.6 


The equation (1.4) can also be written in the form 


7 ] = . - 
d (=) —~ [(E*,) di, + (Ee) dy] + a d (8 *x1X2) 


af +8 BH] 


and determines the adjusted pressure 7, 


Sw dw (1.44) 
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The above results can be stated as 


THEOREM 1. The necessary and sufficient conditions that two axially 
symmetric flows of a viscous liquid corresponding to (%, Xa, Par 2a) 
(a = 1, 2) may be additive are given by the equations (1.5), (1.6) above. 
Equation (1.4) may be used for determining the adjusted pressure 7. 


§2. Self-additive motions. 


From the above considerations it follows that the flow given by (#, x, 
p, 2) is self-additive, i.e., can be added to itself, if and only if, 


8 (ob, & x) 

5(x,o) (2.1) 
and 

4, a a 

Sixx) oe xe PY ~ FRI =O (2.2) 


The equation for determining the adjusted pressure = is easily obtained from 
(1.4) in the form 


7 2 2 l He 24,2 
d (®) = 52 (EW) df + 52 d(x") 


_ 1 Sb\2 (de \2 
alae {(s,) + (x) $1 (2.3) 
From (1.4) we have in the case of self-additive motions, also 

5(F, #) _ 

5(x,a)° (2-4) 
where 

F=" + q’, 

; +q 

showing that 

2 fo” 4. tie 

+ | @*) g () (2.5) 


The equations (2.1) and (2.2) are effectively the same as those obtained 
by S. Chandrasekhar and Kevin H. Prendergast.® 


From (2:1) we get 
Xx = IP (2.6) 
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and on using (2.5) and (2.6) in the equation (2.3) we obtain 
Bb + a SWS WH =8' WH) (2.7) 


THEOREM 2. The class of axially symmetric steady flows with the 
velocity vector 


| ob > 
lx 


2. sae 
a 3a s&o*™ 


is characterized by the conditions (2.6) and (2.7). 


In the case of Beltrami flows with the streamlines coincident with the 
vortex lines, we have F = a constant and g(¥) of (2.5) becomes a constant. 


COROLLARY: The class of axially symmetric, steady, Beltrami flows is 
characterized by the condition 


E*s + f(h) f' (4) = 0 (2.8) 


The equation (2.8) has been completely integrated by S. Chandrasekhar® 
in the case where f(%) =a linear function = af. (An additive constant 


B has no physical consequence.) He has dione that the general solution 
b (x, &) of 


+a?) =0 (2.8 a) 
can be expressed in terms of the fundamental solutions 
on = [AnJnis/2 (ar) + BrJ-nisy2) (ar)] 77*/? Py'*/ (cos 4) 
n= 1, 2, 3, (2.9) 


where r cos ?=-x and rsin?@=@&, An, Bn are arbitrary constants and 
P»°/2 (cos @) denotes the ultraspherical polynomial.* The velocity in such 
a case [f(¥%) =a linear function] of axially symmetric Beltrami flows can 
therefore be represented as any linear combination of the vectors 


( | ddn , cot O ddbn 
r or r2 80 


e 93¢n iI st j ai adn, 


r ér rs rsin 0 


(ne = 1,2, 3, 
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It was shown by P. L. Bhatnagar* that the equation (2-7) can be inte- 
grated on taking g’(¥) =a constant = 8, in addition to f’ (#) =a. 
Then (2:7) becomes 


fe 6 eas or (2.11) 


5x? 


and the solution of this equation is 


a&?-+ solution of (2.8 a). 


2 
< 


The velocity of the flow is given in this case by adding the vector 


mee 8 > 


2 fe to rsin 9 ing 
to any linear combination of the vectors (2.10). 


SUMMARY 


Conditions are noted for the addition of two axially symmetric motions 
of viscous liquids in which the velocities have toroidal as well as poloidal 
components and characterization relations are obtained for the class of 
self-additive axially symmetric motions. 
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